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In order to investigate the orientational ordering of molecular dipoles and the associated electronic 
properties, we studied the adsorption of chlorogallium phthalocyanine molecules (GaClPc, Pc = 
C 3 2N 8 Hi 6 - 2 ) on Cu(lll) using the X-ray standing wave technique, photoelectron spectroscopy, 
and quantum chemical calculations. We find that for sub-monolayer coverages on Cu(lll) the 
majority of GaClPc molecules adsorb in a 'Cl-down' configuration by forming a covalent bond to the 
substrate. For bilayer coverages the XSW data indicate a co-existence of the 'Cl-down' and 'Cl-up' 
configuration on the substrate. The structural details established for both cases and supplementary 
calculations of the adsorbate system allow us to analyze the observed change of the work function. 



PACS numbers: 68.49.Uv, 68.43.Fg, 79.60.Fr 

The adsorption of organic semiconductor molecules has 
been in the focus of numerous experimental and theoret- 
ical investigations - many of them addressing the subtle 
interplay of electronic and structural properties. Early 
studies [1], which show that the energy levels of organic 
semiconductor/metal interfaces can exhibit large devia- 
tions from the Schottky-Mott relation, conveyed the sig- 
nificance of interface dipoles. Until today and despite the 
ubiquity of this concept in the field of organic materials, 
the origin of the interface dipole often remains vague. 

To establish a better understanding of the energy level 
alignment at the interface one should not neglect effects 
related to the molecular structure of organic adsorbates: 
Planar molecules such as Fi 6 CuPc [2], PTCDA [3, 4] 
or pcntaccnc derivatives [5], for example, can distort 
upon adsorption due to the interaction with the sub- 
strate and therefore exhibit an induced molecular dipole. 
Non-planar molecules such as TiOPc [6], SnPc [7-9], 
SubPc [10, 11] and VOPc [12], which may adsorb in 
different orientations, form layers with at least partially 
aligned dipole moments. Hence, for this class of systems 
the orientational order on the surface is a quantity which 
strongly influences the interface dipole. In particular, it 
has been shown that depending on the orientation a layer 
of molecular dipoles p with an area density iVdip can shift 
the vacuum level (VL) in either direction and, therefore, 
increase or decrease the work function <£> of the sample 
according to [13] 

A$ dip = ±epN dip /e e, (1) 

where e is the effective dielectric constant of the mono- 
layer. An experimentally and theoretically challenging 
model system of non-planar organic molecules with a sig- 
nificant dipole moment, for which these effects can be di- 
rectly studied, is chlorogallium phthalocyanine (GaClPc, 



Fig. 1) [14]. 

In this letter, we present a detailed study on the bond- 
ing and orientational ordering of GaClPc on Cu(lll) 
surfaces using the X-ray standing wave (XSW) tech- 
nique [15], ultra-violet photoelectron spectroscopy (UPS) 
and density functional theory (DFT) based calculations. 
While XSW data are taken to determine exact the atomic 
positions along the surface normal and, thereby, also the 
orientation of the molecules, the UPS measurements re- 
veal how the adsorbate affects the electronic energy lev- 
els. In particular, our results demonstrate that the ob- 
served VL shift can be modeled very well if one accounts 
for the adsorption geometry and the orientation of the 
molecules. 

The XSW experiments were carried out at beamline 
ID32 of the European Synchrotron Radiation Source 




FIG. 1: Structure of the free GaClPc (Pc = C 3 2N 8 Hi 6 " 2 ) 
molecule as obtained by geometry optimization using the 
Gaussian03 program package. The results show that the Ga 
(red) and CI atoms (orange) are located above the molecular 
plane with a Ga-Cl bond length of 2.21 A. The Pc group it- 
self is non-planar with the two inequivalent N species (blue) 
0.50 A and 0.58 A below the Ga atom. 
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(ESRF) in Grenoble. Using the (lll)-rcflection of the 
copper substrate we have taken data in back-reflection 
geometry [2, 4], corresponding to the lattice plane spac- 
ing d — 2.08 A with photon energies around -©Bragg = 
2.97 keV. The UPS experiments for different cover- 
ages were performed with a photoelectron spectrometer 
equipped with a standard He I light source (photon en- 
ergy 21.2 eV) in our home laboratory. The Cu(lll) crys- 
tal with a small mosaicity was prepared under ultra-high 
vacuum conditions by repeated cycles of argon sputter- 
ing and annealing (base pressure 4 x 10~ 10 mbar). The 
GaClPc material was purchased from Sigma- Aldrich, pu- 
rified by gradient sublimation, and thoroughly degassed. 
With a quartz microbalance typical deposition rates of 
0.2 A/min were realized. For all experiments reported be- 
low the substrate temperature was carefully monitored. 
To model the electronic structure of the molecules on 
the Cu(lll) surface, we employed slab- type DFT based 
band structure calculations using VASP [1 6]; more details 
on the theoretical methods and the chosen geometry are 
given in the supplementary material [17]. 

The core-level spectra of GaClPc on Cu(lll) show 
several signals that are suitable for XSW experiments 
(Fig. 2). By measuring photoemission spectra related to 
all atomic species in the molecule, namely C(ls), N(ls), 
Ga(2px/ 2 ), and Cl(2s)/C1(KLL) Auger, we achieve a de- 
tailed electronic and structural characterization of the 
system. 
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FIG. 2: Background corrected core- level spectra measured 
with a hemispherical electron energy analyzer for a sub- 
monolayer of GaClPc on Cu(lll). The intensity of the peaks 
relative to the substrate signal has been used to determine 
the surface coverage in monolayer equivalent units (ML). 

The bonding distances and orientation of the molecule 
have been derived from the characteristic variation of the 
photoelectron yield Y p around the Bragg condition. The 
XSW data, which can be analyzed according to a well- 
established procedure using the experimental reflectivity 
R of the sample, yield two structural parameters: the 
coherent position Pjj and coherent fraction fjj, both of 
which can be determined by a least-square fitting rou- 
tine [17]. While Pjj gives the average position of the 
atoms relative to the diffraction planes, fjj yields the 
spread of positions within the ensemble. 

Figure 3 shows XSW results for two different cover- 
ages, which were derived from series of core-level spectra. 



Unlike planar organic molecules [2, 4], the data corre- 
sponding to the different atoms exhibit pronounced dif- 
ferences. For sub-monolayer coverages (0.8 ML, Fig. 3a) 
our results indicate high orientational order within the 
molecular ensemble on the surface. To convert the coher- 
ent positions given in Fig. 3a, i.e. P#(Ga) = 0.03, Pjj(C) 
= 0.14, P ff (N) = 0.27, and P H (CV)= 0.90, into the (av- 
erage) bonding distances djj one must take into account 
the modulo-d ambiguity of the XSW technique [15]. In 
principle, two different configurations can be considered 
for low GaClPc coverages: Either a scenario 'Cl-up' with 
interatomic distances which deviate strongly from the gas 
phase structure (e.g. a hypothetical Ga-Cl bond length of 
only 1.83 ± 0.07 A compared to 2.21 A in the gas phase), 
or a scenario 'Cl-down' with a corresponding Ga-Cl bond 
length of 2.33 ± 0.07 A and a vertical N-Ga distance of 
0.50 A, i.e. values that are close to the interatomic dis- 
tances in the gas phase. DFT calculations for the iso- 
lated molecule, which show that large molecular distor- 
tions of GaClPc are energetically very unfavorable, can 
be used to rule out the 'Cl-up' case. The obvious con- 
clusion that the molecules adsorb - similar to other Pc 
derivatives [10, 11] - predominantly in a 'Cl-down' con- 
figuration implies that the molecular dipoles contribute 
with A$dip < to the overall VL shift. Having estab- 
lished the bonding distances for the low coverage (Fig. 3c) 
with the corresponding error bars (Fig. 3d), we see that 
the Cu-Cl layer spacing of 1.88 A coincides with exper- 
imental results for chlorine adsorption on Cu(lll) [18]. 
This agreement suggests a covalent bonding of the CI 
atoms to the substrate and, moreover, a hollow adsorp- 
tion site as reported in Ref. [18]. 

Complementing this picture, we find increasing mag- 
nitudes of fjj for carbon (0.23), nitrogen (0.41), gallium 
(0.51), and chlorine (0.67), which indicate surprisingly 
different spatial spreads. Therefore, we studied the in- 
fluence of thermally activated librational motions of the 
GaClPc molecules in the 'Cl-down' configuration. Spe- 
cific model calculations based on a statistical description 
of the ensemble (presented in the supplementary mate- 
rial [17]) were performed to refine the analysis. The 
time-averaged spread of positions resulting from our sim- 
ulations gives XSW parameters fjj and Pjj, which - by 
comparison with the experimental values - allow us to 
determine the libration amplitude. The analysis shows 
that the atomic species within the GaClPc molecule are 
affected very differently and that even for small ampli- 
tudes there are strong changes of the coherent fractions 
for C and N. The relative magnitude of the coherent frac- 
tions is consistent with a libration amplitude (rms) of 
17 ± 3° around the essentially fixed Cl-atom. 

For higher GaClPc coverages (1.4 ML, Fig. 3b) the sit- 
uation becomes more complex. The relatively small cov- 
erage dependence found for the CI (and Ga) data con- 
trasts with significant changes for the N (and C) signal. 
In particular the very different values Pr-f(N)=0.49 and 
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FIG. 3: (a),(b) XSW data for GaClPc on Cu(lll) taken for two different coverages in back-reflection geometry at elevated 
temperature. The photon energy has been scanned around i^Bragg = 2.97 keV to measure the photoelectron yield (circles) and 
reflectivity (triangles). Least-square fits to the data (solid lines) give the coherent fraction fn and coherent position Ph that 
are related to the adsorption geometry, (c) Average bonding distances dn for GaClPc as derived from the data in (a), (d) 
X 2 -confidence map for fits to the XSW data in (a). The la-, 2a-, and 3<7-contour levels show the statistical significance of the 
data analysis [4]. 



/#(N)=0.13 for nitrogen cannot be explained with only 
one molecular orientation. Model calculations, which 
take into account that the XSW parameter fn and Pr 
refer to a coherent average over all inequivalent sites in 
the bilayer [15], i.e. 

, N 

f H e 2 " p » = - {f n AfW " /d + e 2 " A7/d ) , (2) 
fe=i 

^down/up d eno ^ m g |- ne position of the k th atom relative 
to the substrate lattice), indicate a co-existence of 'CI- up' 
and 'Cl-down' configurations on the surface. The adsorp- 
tion of GaClPc molecules with opposite orientation in the 
second layer - a phenomenon reported also for other Pc 
derivatives [G, 9] - could be confirmed also by metastable 
atom electron spectroscopy (MAES) measurements [19]. 
Within our coverage series, we verified that a deposi- 
tion of more than 4 ML results in atomically disordered 
systems with incoherent XSW signals {fn = 0) for all 
elements. 

In order to study how the orientational order of the 
molecules affects the VL shift we measured the work 
function of the system for increasing GaClPc coverage 
(Fig. 4a). At first, the results obtained from the pho- 
toemission spectra reveal a rapid decrease below 0.2 ML 
as expected for the 'Cl-down' configuration. The con- 
tinuous, yet slightly retarded VL shift results in a mini- 
mum of A<& = — 0.34eV at the monolayer coverage. To 
overcome possible kinetic limitations we deposited more 
material (resulting in a gradual increase of $ for cover- 
ages of more than 1 ML due to preferential 'Cl-up' nucle- 



ation) and applied a short annealing of the multilayer at 
300 °C. This procedure yields a coverage of slightly less 
than I ML and a shift of A$ = -0.60 eV. Because DFT 
calculations, which we have performed using the adsorp- 
tion geometry shown in Fig. 3c, predict a similar decrease 
of A$ = — 0.55eV for a comparable coverage, we may 
relate the electronic properties of annealed monolayer to 
the low-coverage XSW data. Due to the orientational 
order we can separate the total change A$ into two con- 
tributions with 

A$ = A* dip + A$ bond . (3) 

Thus, we distinguish A$dip, which is related to the 
dipole moment of the molecules according to Eq. (I), 
and A$bond, which contains the effect of the molecule- 
metal interaction. To quantify the dipole contribu- 
tion we have calculated a monolayer of GaClPc with- 
out metal substrate. Again using the adsorption ge- 
ometry shown in Fig. 3c we obtain A$di P = — 0.30eV. 
Because of the large substrate-Pc distance and the cor- 
respondingly small Pauli repulsion the associated value 
A^bond = — 0.25eV essentially reveals the effect of the 
charge rearrangement within the Ga-Cl group and be- 
tween the CI and Cu atoms [17]. 

These conclusion are corroborated by the interpreta- 
tion of the valence band spectra shown in Fig. 4b. The 
data obtained for different coverages show that neither 
before nor after the annealing procedure prominent in- 
terface states - as they were reported for the adsorption 
of other planar [20] and non-planar [21] Pc molecules di- 
rectly below the Fermi energy - can be observed. This 
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FIG. 4: (a) Coverage dependent work function of GaClPc 
on Cu(lll). (b) Valence band region measured by UPS for 
coverages labeled correspondingly in the top panel. In these 
spectra peak A belongs to the Cu3d band, A' is a satellite 
of A due to excitation by parasitic Hel^ radiation, peak B is 
the Shockley surface state of the clean Cu(lll), and peak C 
is the HOMO-derived state of the multilayer film. 



absence of characteristic spectral features in the HOMO- 
LUMO gap [22] illustrates the weak interaction between 
the substrate and the Pc ring in the 'Cl-down' orienta- 
tion. We note that due to the photoemission selection 
rules and the orientation of the molecules no HOMO- 
derived states can be seen in the normal emission spec- 
tra of both monolayer films shown Fig. 4b - in con- 
trast to the disordered multilayer film, which exhibits 
a pronounced HOMO-related peak centered at —1.42 eV. 
Moreover, the spectra illustrate the non-ideal growth of 
GaClPc at 80 °C and the significant re-ordering caused 
by the annealing procedure. The strong quenching of 
the Shockley surface state, which we find for the annealed 
film, clearly indicates a more uniform substrate coverage. 

In summary, the XSW data together with the UPS 
measurements and the DFT calculations give a conclusive 
picture of the molecules' orientation and bonding. For 
sub-monolayer coverages at sufficiently high substrate 
temperature the vast majority of GaClPc molecules ad- 
sorb in the 'Cl-down' configuration with a Cl-Cu layer 
spacing of 1.88 A that indicates covalent bonding of the 
chlorine atoms. The interaction with the substrate gives 
rise to subtle deviations from the gas phase structure: 
In agreement with the DFT calculations [17], the Ga-Cl 



bond length of 2.33 A measured by XSW is slightly larger 
than in the gas phase, while the average position of the 
carbon atoms in the experiment and the calculations in- 
dicates a certain bending of the molecule [17]. Moreover, 
the coherent fractions observed in the experiment can be 
related to thermal librations of the molecules. Due to 
the orientational order the molecular dipoles contribute 
significantly to the vacuum level shift and the result- 
ing interface dipole. Detailed calculations for the an- 
nealed monolayer support this conclusion quantitatively 
with A$dip = — 0.30eV. The influence of intermolecular 
interactions is demonstrated by the coexistence of 'Cl- 
down' and 'Cl-up' configurations found for deposition at 
80 °C. Interestingly, a subsequent annealing allows a con- 
trolled rc-ordcring of the molecular dipoles and a change 
of the corresponding electronic properties. Furthermore, 
we believe that the adsorption behavior of GaClPc and 
its impact on the electronic structure represent impor- 
tant characteristics of non-planar organic molecules. 
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molecular orbital 



[2 
[•3 

[5 

[6 
[7] 

[8] 

[0 

[10 
[11 

[12 
[13 

[14 
[15 
[10 

[17 
[18 
[19 
[20 

[21 

[22 



